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We calculate the size distribution of quark nuggets, which could be formed due to 
first order QCD phase transition in the early universe. We find that there are a large 
number of stable Quark Nuggets which could be a viable candidate for cosmological 
dark matter. 

1 Introduction 

A first order quark-hadron phase transition in the early universe at a critical tem- 
perature T c ~ 100 — 200 MeV, lead to the formation of quark nuggets (QN), made of 
u, d and s quarks Under certain circumstances the primordial QN's will survive 
till the present epoch. The central theme of this work is the candidature of these 
quark nuggets as the baryonic component of the dark matter. This possibility leaves 
the results of big bang nucleosynthesis unaffected and does not invoke any exotic 
physics 

One of the significant issues in this context is the stability of these primordial 
QN's on a cosmological time scale. This question was first addressed by Alcock and 
Farhi |J who argued that due to baryon evaporation from the surface, a QN, even 
with the largest possible baryon number, will not survive till the present epoch. 
Madsen et al || then pointed out that as the evaporation makes the surface of 
the nugget deficient in u and d quarks, further evaporation is suppressed. They 
came to a conclusion that QNs with initial baryon number Nb > 10 46 could well 
be stable against baryon evaporation. Later Bhattacharjee et al [[F] found, using the 
chromoelectric flux tube model, that QN's with baryon number larger than 10 42 , 
would survive against baryon evaporation. 

In spite of these efforts, not much emphasis has been put towards the study of 
the size-distribution of the QNs. The size-distribution of the QNs is very important 



in the context of their candidature as dark matter as it tells us the most probable 
size of a QN. The calculation of the lower cut-off in size tells us the minimum size 
and the baryon number content of a QN that we should look for. We will carry out 
these studies in the cosmological QCD phase transition scenario. 



2 The size-distribution of quark nuggets 

The evolution of the cosmological scale factor during the quark-hadron phase tran- 
sition epoch is given by, 
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where r = g q / gh, t c = ^m^jSirB, ti is the time when phase transition starts and 
B is the bag constant. (For details and explanation of all the terms see ref. 0). 

In the coexisting phase, the temperature of the universe remains constant at T c . 
In the usual picture of bubble nucleation in a first order phase transition scenario 
hadronic matter starts appearing as individual bubbles. With the progress of time, 
more and more hadronic bubbles form, coalesce and eventually percolate to form an 
infinite network of hadronic matter which traps the quark matter phase into finite 
domains. The time when the percolation takes place is usually referred to as the 
percolation time t p , determined by a critical volume fraction f c , (f c = f(t p )) of the 
quark phase. 

In an ideal first order phase transition, the fraction of the high temperature 
phase decreases from the critical value f c , as these domains shrink. For the QCD 
phase transition, however, these domains should become QN's and as such, we may 
assume that the lifetime of the mixed phase tf ~ t p . The probability of finding such 
a domain of trapped quark matter of co-ordinate radius X at time t p with nucleation 
rate I(t) is given by [[|, 
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where X(t p ;t) is the coordinate radius of a bubble, at time t p , which nucleated at 
time t. 

For convenience, we define a new set of variables z = XR(ti)/vt c , x = t/t c and 
r(x) = R(x) / R(xi); v is the radial growth velocity of a bubble. Then 
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In order to find the minimum size and the size-distribution of the QNs we follow 
the procedure of Kodama et al 0. The distribution function, in terms of z, is given 
by!,! 
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The solution of the equation F(a) = gives the minimum size, a, of the quark- 
nugget. The number of nuggets per unit volume is then 

roo poo f}3(f.\ 

n Q = R- 3 (t p ) / F(X)dX = R~ 3 (t p ) / ±£g-f(z)dz (7) 
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The volume of each quark nugget is given by ^7i(zvt c ) 3 . Since visible baryonic 
matter constitutes only ten per cent of the closure density (fig — 0.1 from standard 
big bang nucleosynthesis), a total of 10 50 baryons will close the universe baryonically 
at T = 100 MeV. We emphasize at this point that these QNs would not disturb the 
standard primordial nucleosynthesis results. Therefore, if we assume that the total 
baryon content of the dark matter is carried by the quark nuggets then, 

N B = 10 50 (100/T(MeV)) 3 = V H ^^-p ^ f(z)z 3 dz (8) 

where Vh is the horizon volume and p is the baryon density inside each nugget. 
We have taken p = 0.15/m~ 3 and v = 0.5 in the present calculation. The above 
equations are then solved self-consistently to obtain a and t p . These values are 
then used to study the size-distribution of the quark nuggets. To calculate the size 
distribution of QNs we have used the nucleation rates proposed by Csernai and 
Kapusta [|1(J. 

Once the baryon density inside the nuggets is known one can easily translate 
from z to ub (baryon number of a particular QN). In fig. [I| we have plotted the 
distribution of QN, f(n B ), as a function of n B . We see that for T c = lOOMeV^ there 
is no quark nugget below n B = 10 46 and above n B = 10 47,5 . For T c = 150MeV there 
is no quark nugget below n B = 10 41 ' 5 and above n B = 10 43 5 . Earlier studies 0] 
have shown that the nuggets having baryon number less than 10 42 will not survive 
till the present epoch. So all the nuggets for T c = lOOMeV will survive and some 
of the nuggets for T c = 150MeV will survive. 




Figure 1: Distribution of QN, /(%), as a function of ub using nucleation rate 
proposed by Csernai and Kapusta. The value of surface tension is 50MeV fmr 2 . 



3 Conclusion 



In this work we have estimated the abundance of quark nuggets in various nucleation 
scenarios with different values of critical temperature. We have found that within 
a reasonable set of parameters QNs may be a possible candidate for cosmological 
dark matter. 

One of us (JA) is grateful to Japan Society for Promotion of Science (JSPS) for 
financial support. 
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